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I. INTRODUCTION 

We have developed over the last years a general 
strategy to calculate physical observables of com- 
pound systems in a nonperturbative framework [THB] . 
This strategy is based on light-front dynamics (LFD), 
a relativistic Hamiltonian formalism advocated first 
by Dirac in 1949 [7 . In its original formulation, the 
state vector of any compound system evolves in the 
light-front time r = t + z, instead of the usual time t. 
Since the physical vacuum in LFD is identical to the 
free vacuum (see e.g. Ref. [8] and references therein), 
it is then natural to decompose the state vector of 
a compound system in Fock components, since this 
decomposition does not include any vacuum fluctua- 
tions, but contains physical (asymptotic) states only. 
The problem of finding the state vector can thus be 
formulated as a TV-body problem. 

For obvious practical reasons, the Fock decomposi- 
tion should be truncated to a finite number of states 
(Fock sectors) . This truncation should be strictly con- 
troled in order to be able to make reliable predictions 
for physical observables, order by order in the Fock ex- 
pansion. This is made possible due to two important 
breakthroughs: 

• The formulation of LFD in a covariant way , 
which enables the strict control of any violation 
of rotational invariance, when the Fock space 
is truncated. This formulation, called covariant 
light-front dynamics (CLFD), has proven to be 
very powerful in the description of relativistic 
properties of few-body systems [ID] . 



• The development of an appropriate renormal- 
ization procedure — the so-called Fock sector 
dependent renormalization (FSDR) scheme 
which enables to calculate regularization scale 
invariant observables order by order in the Fock 
expansion g]. 

In the simplest, two-body, Fock space truncation, 
our formalism is equivalent to summing the irreducible 
block — the fermion self-energy calculated in the sec- 
ond order of perturbation theory — to all orders in 
the chain approximation. This equivalence is caused 
by the fact that all the chain type contributions are 
restricted to the two-body Fock sector. Such a result 
is a direct consequence of our FSDR scheme and the 
corresponding renormalization conditions. 

The first nontrivial calculation corresponds to the 
three-body Fock space truncation which incorporates, 
in the Yukawa model for instance, fluctuations of the 
state vector involving one fermion (/), one fermion 
and one boson (fb), one fermion and two bosons (fbb) 
Fock sectors. This calculation includes overlapping 
type (divergent) diagrams summed to all orders in the 
coupling constant. 

Within the FSDR framework, the first calculation 
of a physical observable — the anomalous magnetic 
moment (AMM) of a fermion in the Yukawa model — 
has been done in Ref. [5], using the Pauli-Villars (PV) 
regularization scheme, as proposed in Ref. [IT]. The 
calculation has shown nice convergence of the results 
as a function of the regularization scale (the PV bo- 
son mass in our calculation) for values of the coupling 
constant a = g 2 /^ of order of 0.2. For stronger cou- 
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pling, a ~ 0.5, some dependence (though rather weak) 
of the AMM on the PV boson mass was detected. 
While this range of the coupling constant values is 
not particularly small (as compared for instance to 
the electromagnetic coupling constant), it shows how- 
ever that the truncation of the Fock expansion was 
not completely under control. 

We detail in the present study an extension of our 
previous approach 0, in order to control, order by 
order in the Fock expansion, the regularization scale 
invariance of physical observables. Our derivation is 
based on the full account of the renormalization con- 
ditions, using the FSDR scheme. We shall see that 
in the truncated Fock space, the bare coupling con- 
stant and counterterms are no more true constants, 
but become naturally dependent on one of the kine- 
matical variables, as already emphasized in Ref. [T^] . 
This dependence is determined unambiguously by the 
renormalization conditions and allows to restore ro- 
tational invariance broken by the truncation, as well 
as the independence of observables on the masses of 
the PV particles, when the latter ones are much larger 
than the physical masses. In Ref. [5] the fact that the 
bare coupling constant and counterterms are, a priori, 
functions of kinematical variables was not taken into 
account in full measure. The calculations were done 
for fixed values of the kinematical variables. 

The kinematical dependence of the bare coupling 
constant and counterterms is intimately linked to the 
Fock space truncation. Its explicit form is strongly 
affected by the Fock space "contents". Thus, in the 
Yukawa model considered in the three-body (/ + /& + 
fbb) approximation, this dependence is quite sizeable. 
It is substantially reduced, when antifermion degrees 
of freedom (d.oi.), namely, the additional three-body 
Fock sector ///, are taken into account. In leading or- 
der of perturbation theory, we show that the inclusion 
of the antifermion d.o.f. results in the independence 
of the bare coupling constant and counterterms on the 
kinematical variables. 

The plan of the article is as follows. We recall in 
Sec. |TT| the general properties of our formalism. We 
discuss in Sec. HIT] the renormalization conditions. In 



Sec. [TV] we obtain a system of renormalized equations 
for the Fock components in the Yukawa model within 
the / + fb + fbb Fock space truncation and calculate 
the fermion electromagnetic form factors. In SecJVjwe 
extend the Fock space by the inclusion of the /// Fock 
sector and discuss the role of antiparticle d.o.f. We 
present our conclusions in Sec. [Vl] The contribution 
of antifermion d.o.f. to the equations for the Fock 
components is given in Appendix |A| 



II. GENERAL FRAMEWORK 
A. Covariant formulation of light-front dynamics 

In the traditional form of LFD, the state vector of 
a compound system is defined on the light-front plane 
t + z = (with c = 1) rather than on the equal-time 
plane t = 0. In order to recover explicitly rotational 
invariance, the state vector is defined, in CLFD, on 
the light- front plane of general orientation uj-x = 0, 
where uj is an arbitrary four-vector restricted by the 
condition ou 2 = [9 , 10 . The traditional form of LFD 
is recovered by using u = (1,0,0, —1). 

The state vector 4>(p) of a particle with the mass 
M should satisfy the Poincare group equations, and 
among them 



p 2 0(p) = m 2 <Kp). 



(I) 



The momentum operator P is decomposed, on the 
light front, in terms of its free and interaction parts: 



(2) 



where, in terms of the interaction Hamiltonian 

H int (x), 



pint 
P 



H int (x)5{u-x) d A x. 



(3) 



According to the general properties of LFD, we de- 
compose the state vector of a physical system in Fock 
sectors. We have schematically 

oo „ 

</>(p) = ^2 / db n <f> n (ki, . . . , k n ;p) 

n=l J 

x (5( 4) (/ci + ... + /c n -p-wT n )|n), (4) 

where \n) is the state containing n free particles with 
the four- momenta k±,...,k n and n 's are relativis- 
tic n-body wave functions, the so-called Fock compo- 
nents. The phase space volume element is represented 
by dD n . All the four- momenta are on the correspond- 
ing mass shells: kf =m 2 , p 2 = M 2 , (ujr n ) 2 = 0. Note 
the peculiar overall four-momentum conservation law 
given by the (5-function. It follows from the gen- 
eral transformation properties of the light-front plane 
lu-x = under four-dimensional translations [TO]. The 
scalar quantity r n is a measure of how far the n-body 
system is off the energy shell (on the energy shell 
r n = 0). It is completely determined by this conser- 
vation law and the on-mass-shell conditions for each 
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individual particle momentum. We get 
2uj-p r n = (s n - M 2 ), with s n = (fci + . . . 4- /c n ) 2 - (5) 



The state In) can be written as 



\n)=rf(k 1 )rf(k 2 )...rf(k n )\Q), 



(6) 



where d) is a generic notation for the fermion and 
boson creation operators. To completely determine 
the state vector, we normalize it according to 



0(p') t 0(p) = 2p o <5 (3) (p'-p)- 



(7) 



With the decomposition Q, the normalization condi- 
tion ([7]) writes 



(8) 



n=l 



where I n is the contribution of the n-body Fock sector 
to the norm. For the particular case of the Yukawa 
model, an explicit formula for I n can be found in 
Ref. [J. 

We shall concentrate, in the following, on systems 
composed of a spin- 1/2 fermion coupled to scalar 
bosons. It is convenient to introduce, instead of the 
wave functions (f) n , the vertex functions T n (which we 
will also refer to as Fock components), defined by 



u{k\)T n u{p) = (s n - M 2 )(f) n = 2u-p r, 



(9) 



where k\ is the four-momentum of the constituent 
fermion. When the Fock space is truncated to order TV 
[i. e. in the sums over n in Eqs. Q and ^ the terms 
with n < N only are retained], it is necessary to keep 
track of the order of truncation in the calculation of 
the vertex function. For this purpose, we will denote 
the latter by Tn\ but omit the superscript (N) when 
it is not necessary. is represented graphically by 
the diagram shown in Fig. [T] 



(n — 1) bosons 




It is convenient to introduce the notation 



Q(p) = 2(u-p)r(f)(p), 



(10) 



where f is the operator which, acting on a given com- 
ponent (p n of 0(p), gives r n (j) n . Q(p) has the Fock 
decomposition which is obtained from Eq. Q by the 
replacement of the wave functions cf) n by the vertex 
functions T n . We can thus cast the eigenstate equa- 
tion in the form [1^ 



(11) 



where H int is the interaction Hamiltonian in momen- 
tum space. This equation is quite general and equiv- 
alent to the eigenstate equation ([!]). It is nonpertur- 
bative. 

The graph technique rules derived in Ref. [10] for 
the calculation of ^-matrix elements in CLFD re- 
fer precisely to the calculation of matrix elements of 
—H mt /(2iTT). A system of coupled eigenstate equa- 
tions for the Fock components of the state vector can 



thus be obtained from Eq. (11) by substituting there 
[via Eq. ( 10)] the Fock decomposition Q and calculat- 
ing the matrix elements of the operator — H int /{2nr) 
in the Fock space. 

Since our formalism is explicitly covariant, the spin 
structure of the wave function cj) n is very simple. In- 
deed, its construction is of purely kinematical na- 
ture [ID] . It should incorporate however o;-dependent 
components. The structure of the two-body compo- 
nents for QED and the Yukawa model was detailed in 
Ref. 0, while the three-body (fbb) component in the 
Yukawa model was constructed in Ref. [5 . For the 
purpose of the present study, we recall here the spin 
structure of the two-body component in the Yukawa 
model: 



u{ki)T 2 u(p) = u(ki 



h 4- M * h 

01 H 02 

uj-p 



u(p). (12) 



The coefficients bi and b 2 are scalar functions deter- 
mined by dynamics. 



B. Fock sector dependent renormalization 
scheme 



FIG. 1: n-body vertex function for the Fock space trun- 
cation of order N, for a physical fermion (double straight 
line) made of a constituent fermion (single straight line) 
coupled to (n — 1) bosons (wavy lines). 



In standard renormalization theory, the bare pa- 
rameters (the whole set of bare coupling constants 
and counterterms) are determined by relating them 
to physically observable quantities. To perform this 
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strategy in practice, two important questions should 
be clarified, when using CLFD. 

(i) The explicit form of the relationship between the 
bare and physical parameters depends on the approx- 
imation used. For instance, in perturbation theory 
the bare parameters are polynomial functions of the 
physical coupling constant. The term with the max- 
imal degree of the coupling constant is determined 
within a given order of perturbation, while those of 
lower degrees are taken unchanged from lower orders. 
This guarantees that calculated physical observables 
are regularization scale invariant in any order of the 
perturbative expansion. In our nonperturbative ap- 
proach based on the truncated Fock decomposition, 
an analogous requirement implies that the bare pa- 
rameters should depend a priori on the Fock sector in 
which they are considered [13]. 

(ii) In order to express the bare parameters through 
the physical ones, and vice versa, one should be able 
to calculate physical observables. In LFD, these ones 
can not depend on the choice of the orientation of 
the light-front plane. Such a situation indeed takes 
place, for instance, order by order in perturbation 
theory, provided the regularization of divergences in 
LFD amplitudes is done in a rotationally invariant 
way 0. In nonperturbative LFD calculations, which 
always imply Fock space truncation, the dependence 
on the light-front plane orientation may survive even 
in calculated physical amplitudes. The identification 
of such amplitudes with observable quantities becomes 
therefore ambiguous. 

The use of our FSDR scheme in CLFD allows us 
to answer both questions. In this scheme, each of the 
original bare parameters has an additional index de- 
pending on the number of particles in the Fock sector 
in which this bare parameter appears. In the Yukawa 
model the fermion mass counterterm, Sm 1 and the 
bare coupling constant go, are thus extended each to 
a whole sequence: 

9o -> 9oi , (13a) 
5m -)• 5mi, (13b) 

with 2 = 1,2,... iV. By definition, g i = and 
5m\ = 0. For I > 2, the quantities goi and 5mi are 
calculated by solving the systems of equations for the 
vertex functions in the N = 2, TV = 3, ... approxima- 
tions successively. 

Besides that, as we shall see in the next section, 
new counterterms which depend explicitly on the ori- 
entation of the light-front plane (i. e. on uj) should 
be introduced in order to restore, if necessary, the ro- 



tational invariance broken by the truncation. In this 
case, such counterterms are also mandatory in order 
to fulfil the renormalization condition. Similarly to 
the "traditional" counterterms, they have Fock sector 
dependence. The full set of rules for the calculation 
of the bare parameters can be found in Refs. [4j [6] . 

We emphasize that the FSDR scheme is a general 
method to make nonperturbative calculations in trun- 
cated Fock space. It can be easily applied to any phys- 
ical system admitting particle counting. The Yukawa 
model studied in the present work has been chosen as 
an illustration of the capabilities of our approach. 

III. RENORMALIZATION CONDITIONS 

Once the bare coupling constant and the mass coun- 
terterms have been identified, one should fix them 
from a set of renormalization conditions. In pertur- 
bation theory, there are three types of quantities to 
determine: the mass counterterms and the bare cou- 
pling constant just mentioned, and the field strength 
renormalization constants. In the on-mass-shell renor- 
malization scheme, the following conditions are used. 
The mass counterterms are fixed from the requirement 
that the two-point Green's functions have a pole at the 
physical masses. The field normalization constants are 
fixed from the condition that the residues of the two- 
point Green's functions at these poles equal 1. The 
bare coupling constant is determined by requiring that 
the on-mass-shell three-point Green's function is given 
by the product of the physical coupling constant and 
the elementary vertex. 

The renormalization conditions in LFD are of 
slightly different form, although they rely on the same 
grounds. The mass counterterm for each physical 
state is fixed from the eigenstate equation ([!]) by de- 
manding that the mass, M, of the physical bound 
state be identical to the constituent mass m. The bare 
coupling constant is determined by relating the on- 
energy-shell two-body vertex function T 2 to the phys- 
ical coupling constant g. Finally, the normalization of 
the state vector is fixed from the condition ([7]). 

In order to set up the relationship between T 2 and 
the physical coupling constant, one needs to discuss 
carefully the renormalization factors of the external 
legs of the two-body vertex function [5] 03] . These 
renormalization factors do also depend on the order of 
truncation of the Fock space. In the Yukawa model, 
this relationship reads 

r< W )( S2 = M 2 ) = g^p ZI > x ^. (14) 
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Eq. ( 14 ) can be interpreted in simple physical terms. 



Each leg of the two-body vertex function contributes 
for a different factor \J~Z to the physical coupling 
constant, where Z is the field strength normaliza- 
tion factor. The initial fermion state is the phys- 
ical state normalized to 1, so that Z = 1 in that 
case. The final boson line should be renormalized by 
a factor \/Zb. In the approximation where fermion- 
antifermion loop contributions are not considered (the 
so-called quenched approximation), we have Z& = 1. 
Finally, the field strength normalization factor of the 
constituent fermion is just the weight of the one-body 
component in the norm of the physical state [5 , i.e. 
Z = Ji in that case, according to Eq. Q. Following 
our FSDR scheme, the normalization factor of the fi- 
nal constituent fermion should correspond to the trun- 
cation of order N — 1 of the Fock space, since there is, 
by construction of the two-body vertex function, one 
extra boson in flight in the final state. 

Under the PV regularization, PV particles are con- 
sidered, in the interaction Hamiltonian, on equal 
grounds with the physical ones. From here it follows 
that each constituent particle line in the two-body ver- 
tex may correspond to either a physical or a PV parti- 
cle. Observable amplitudes are described by diagrams 
with physical external legs only. For this reason, the 
renormalization condition (14) should be imposed on 



the two-body vertex function with constituent lines 
corresponding to the physical fermion and boson. 

The condition ( 14 ) has an important consequence: 
the two-body vertex function at S2 = M 2 should be in- 
dependent of the orientation uj of the light-front plane. 



With the spin decomposition (12), this implies that 
the component b 2 at S2 = M 2 should be identically 



zero: 



bi N) (s 2 = M^) = 0. 



(15) 



If Eq. (15) is satisfied, Eq. (14), in the quenched ap- 



proximation, turns into 



5 2 = M ) EE 



(N-l) 



(16) 



While the property Jl5| ) is automatically verified in 
the case of the two-body Fock space truncation, pro- 
vided one uses a rotationally invariant regularization 
scheme [3 [4] , this is not guaranteed for higher order 
truncations. Indeed, nothing prevents T 2 to be u- 
dependent, since it is an off-shell amplitude, but this 
dependence must completely disappear on the energy 
shell. It would be so if no Fock space truncation has 
been done. The latter results in some ^-dependence of 



T2 even on the energy shell. This immediately makes 
the general renormalization condition (14) ambiguous, 



since its right-hand side is a;- independent. 

Another consequence of the truncation of the Fock 
space is the fact that the components 61,2(52 = M 2 ) 
are not constants. Indeed, depend a priori on 
two kinematical variables. For practical purposes, we 
can take the usual longitudinal momentum fraction x 
and the transverse (with respect to the three-vector 
u>) momentum Ri. They are defined by 



uj-ko 



x = 



uj-p 

R ± = k 2± - a?p ± , 



(17a) 
(17b) 



where k 2 refers to the momentum of the boson in the 
two-body Fock sector. Note that R ± = — (k 2 — xp) 2 
is a relativistic invariant. We have therefore b\ j2 = 
b 1 , 2 (R ± ,x). 

We denote by m (fi) the constituent fermion (boson) 
mass. The on-shell condition 



R 2 +ra 2 R 2 

S2 = + 

1 — x : 



= M 2 



(18) 



can be used to fix one of the two variables, say in 
the non-physical domain (for M = m): 



R ± = R* ± (x) ee i^/x 2 m 2 + (l-x)/i 2 , 



(19) 



so that b[ N 2 \s 2 = M 2 ) ee (R* ± (x) , x) calculated in 
the truncated Fock space depend on x [an example of 

a particular form of the function b^jR* (x),x) for 
N = 3 is given in Sec. |V| Eq. (57)], whereas the 
conditions ( [15] ) and ( p~6) should be valid identically, 
i.e. for any value of x. 

In order to enforce the condition ( fl5| ), we should 
introduce an appropriate counterterm which depends 
explicitly on the four- vector w 0. It originates from 
the following additional structure in the interaction 
Hamiltonian: 



8HT = -Z u i/>- 



jr 171 <fi 



iuj-d 



(20) 



where Z u is just the new counterterm, ijj'((p') is the 
fermion (scalar boson) field, being a sum of the corre- 
sponding physical and PV components, and l/{iu-d) 
is the reversal derivative operator. In the Yukawa 
model within the three-body approximation the con- 
tribution p0| ) is enough to make all renormalization 
conditions self-consistent. 



6 



In the truncated Fock space, according to the FSDR 
rules, splits into a sequence of Fock sector depen- 
dent contributions Z<y , analogously to the other bare 



x-dependence of the bare parameters should become 
weaker. We will see an indication of that behavior in 



parameters [see Eqs. ( 13 )]. For the truncation of order 

N, it is supposed that all Z^'s with I = 1, 2, . . . N- 1 
have been already known from lower order trunca- 
tions, so that we have to determine the "senior" coun- 
terterm only. The enforcement of the condi- 



tion (15), for any x, by an appropriate choice of 



the counterterm Z^ implies that the latter should 
a priori depend on x, i.e. Z^ = Z^ N \x). If no 
Fock space truncation occured, we would get the ex- 
act equality Z w = const = 0, like e.g. in perturba- 
tion theory. The same happens for lowest order Fock 
space truncations because of their triviality. Thus, 
zffl = by definition. Then, in the two-body ap- 

(2) 

proximation, Z^ ' is also zero, provided the PV regu- 
larization is used [4] . Nonzero and ^-dependent coun- 
terterms Z^ N \x) appear, starting from N = 3. 

Following the above discussion, the enforcement of 
the condition ( 16 ) induces also a unique dependence of 



9on = 9qn{x) as a function of the kinematical variable 
x. 

The fact that, in order to satisfy the renormaliza- 
tion conditions, the bare parameters must depend on 
the kinematical variable x, is crucial to obtain results 
which are finite after the renormalization procedure in 
the truncated Fock space is applied. In Sec. |IV C[ the 
stability of our results relative to the value of the reg- 
ularization scale, if the latter reasonably exceeds the 
physical masses, will be confirmed numerically with 
high precision. 

At first glance, the x-dependence of the bare param- 
eters seems, at least, unusual. However, it is a natu- 
ral consequence of the truncation. Of course, the bare 
parameters in the fundamental non-truncated Hamil- 
tonian are true constants. After truncation, the initial 
Hamiltonian is replaced by a finite matrix which acts 
now in a finite Fock space. But it turns out that the 
modification of the Hamiltonian is not restricted to a 
simple truncation. Indeed, to preserve the renormal- 
ization conditions, the bare parameters in this finite 
matrix become ^-dependent. This x-dependence can- 
not be derived from the initial Lagrangian. It appears 
only after the Fock space truncation. 

Our truncated Hamiltonian with x-dependent bare 
parameters is a self-consistent approximation to the 
initial fundamental Hamiltonian. One expects that 
the approximation becomes better, when the number 
of Fock components increases. At the same time, the 



Sec. |VB[ We emphasize that there is no any ambi- 
guity in finding the bare parameters, in spite of their 
x-dependence. They are completely fixed from the 
renormalization conditions. 



IV. YUKAWA MODEL IN THE f + fb+ fbb 
APPROXIMATION 



We apply our general strategy to calculate some 
physical observables for the Yukawa model in the trun- 
cated Fock space including sectors with one single 
fermion, one fermion plus one boson, and one fermion 
plus two bosons. Previously, we considered the same 
physical system to calculate the fermion AMM [5 J , but 
without x-dependent bare parameters. The AMM (as 
well as any calculated observable) depends on the reg- 
ularization parameters which are the two PV fermion 
and boson masses mi and /ii, respectively. In case 
of a proper renormalization scheme it must tend to a 
fixed finite value, when both PV masses become much 
greater than the characteristic physical mass scale. 
In Ref. [5] we first took the limit mi — >• oo analyti- 
cally, just on the level of the equations for the Fock 
components, and then studied the dependence of the 
AMM on the remaining PV mass \i\ numerically. We 
found, at relatively small values of the coupling con- 
stant (a ~ 0.2), rather good numerical stability of the 
AMM as a function of fii, when ji\ ^> m, fi. 

At larger coupling constants (a ~ 0.5), we observed 
weak but sizable growth of the AMM with the in- 
crease of fi\. A possible reason of this uncanceled 
/ii-dependence of the AMM is the fact that we used 
constant, i. e. x-independent, bare parameters. We 
shall show in this study that taking into account x- 
dependence of the bare parameters in the truncated 
Fock space allows to remove completely any depen- 
dence of observables on the regularization parameters, 
even for rather large values of the coupling constant. 
We shall thus calculate not only the AMM, but both 
electromagnetic form factors as a function of the mo- 
mentum transfer squared. 



A. Equations for the Fock components 

As shown in Ref. [5j, the system of equations for 
the three vertex functions can be reduced to a closed 
matrix equation which involves the two-body vertex 
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function r 2 only. This equation is shown schemati- 
cally in Fig. [2] Each factor at the vertices is taken 
according to the FSDR scheme prescriptions. The 
factors calculated in the three-body approximation, 
namely, #03 and appear in the amplitudes involv- 
ing the one-body vertex function Ti only. Each of the 
other boson emission and absorbtion vertices in Fig. [2] 
brings the factor g 2, since there exists one boson in 

flight at the time moment corresponding to the ver- 

(2) 

tex and Z^ ' =0. The mass counterterm contribution 
(61712) denoted by the cross appears, in the equation 
for T 2 , within the two-body sector only. Analogous 
contributions inside the three-body sector are absent, 
because they correspond to the factor 5mi which is 



zero. 



r 3 



<)():>, 



u ui.p 




sort of constituent. Each index may take two values, 
and 1, corresponding to a physical or a PV particle, 
respectively. The substitution T 2 — > T 2 J means the 
substitutions b\ )2 — > b\ 3 2 for each of its spin compo- 



nents defined by Eq. (12). The particle masses are 



denoted by and with mo = m and /i = /i, by 
definition. 



We shall consider in this study two different orders 
in taking the infinite limit for the masses of the PV 
particles. One of them, already considered in Ref. pj, 
corresponds to the case, when the mass of the PV 
fermion, mi, is first set to infinity (analytically), while 
the PV boson mass, /ii, is kept finite. In the present 
section, we re-formulate the system of equations for 
the Fock components in this limit, taking into account 
the ^-dependence of the bare parameters. We then 
study numerically the dependence of our results on 
/xi, when the latter is much greater than the physi- 
cal masses. In addition, we discuss in Sec. |IVD| the 
opposite order of limits, when [i\ is first set to infinity. 



FIG. 2: Equation for the two-body Fock component in the 
/ + fb + fbb approximation. 



In order to take the limit mi — >• 00, it is convenient 
to introduce a set of functions, h\ and H%, which re- 
main finite in this limit. They are given by [5] 



The fermion-boson loop (see the third diagram on 
the right-hand side of the graphical equation in Fig.[2| 
is nothing else than the two-body fermion self-energy. 
In CLFD, it can be represented as the following gen- 
eral decomposition [3) [5]: 



S(» = 9o 2 2 



m uj-k 



, (21) 



where A, £>, and C are scalar functions given in 
Ref. [5 . Note that under the PV regularization 
scheme we have C(k 2 ) = 0. 

Since each of the two constituent lines of the two- 
body vertex may correspond to either a physical or 
a PV particle, we have to distinguish four types of 
T 2 , depending on its "contents": (i) physical fermion 
and physical boson; (ii) physical fermion and PV bo- 
son; (Hi) PV fermion and physical boson; (iv) two PV 
particles. It is convenient to supply T 2 by the two in- 
dices, i (for fermion) and j (for boson), reflecting the 



b\ J = 

m 



*2 = 



m , 



Hi-ii-x+mrt 



in 



2(1 -x) 



(22) 



In terms of the functions hj and Hf, the equation 
shown in Fig. [2] corresponds to the following system 
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of integral equations [5] 1 : 



and re-define the bare parameters by 



hi(R ± ,x) = g'^+g' 2 Kihi(R ±1 x)+Kihi(R ±1 x) 



+g' 2 i 3 (R±,x), 



(23a) 



h{{R ± ,x) = g' 2 [-K J 3 hi(R ± ,x) + K{h{(R^x) 

+g' 2 i{{R u x\ (23b) 
H'(R ± ,x) = g' m [(2-x) + Z' UJ {l-x)] 

+g' 2 [K(Hi(R ± ,x) + K 3 2 Hi(R ± ,xj 
+g ,2 Ij>(R ± ,x), (23c) 



H{(R ±1 x) = 



003 



+g' 2 -KiH 3 (R ± ,x) + KiHi(R ± ,x) 



+g' 2 I{(R^x). 



(23d) 



where 



with 



003 



= 003^0 



i + g 2 z 



z' = 



2Z 



(3) 



003 



Ql 

Q!o' 



7 (2)_ B(m 2 ) 

^2 — Z 0, 

m 



(26a) 
(26b) 



(27) 



In the above equations we use the following notations. 
The coefficients K{_ 4 are defined by 



z = 2m [A!(m 2 ) + B'{m 2 )]. 



(28) 



K{ = 



m 



— si 



R j = A r (Sl) +_gr(gl) 



m[A(si) + Br(«i)] 



m z — si 



Br(si) 



where the subtracted fermion self-energy functions are 
A r (s!) = A(s 1 )-A{m 2 ), Br(si) = B(s 1 )-B(m 2 ) 1 and 
their argument 



jl i /-i \ 2 ^" ^ 2 

si = h (1 — x)m — // • 



3 ' 



(24) 



For convenience, we introduce the quantity 



9 ' 2 = 



i + g 2 z 



(25) 



The primes at A and B denote their derivatives 
over s\. The quantities oto and q>\ are, respectively, 
the physical and PV components of the one-body 
(fermion) vertex, introduced in Ref. [5]. We do not 
need to calculate them explicitly in our study. We 
can therefore include them into the definitions of g' 03 
and Z' u . The quantity I 2 , multiplied by g 2 , is noth- 
ing else than the norm of the two-body sector, cal- 
culated for the two-body (/ + fb) Fock space trunca- 
tion. Transforming the equation for T 2 to the system 
of equations ( 23 ) , we took into account the values of 



the bare coupling constant #02 and the mass countert- 
erm 8m2, obtained from the previous calculations g] 
within the / + fb truncated Fock space 2 : 



002 = 



l-^ 2) ' 
5m 2 = g 2 2 [A(m 2 )+B(m 2 )]. 



(29a) 
(29b) 



The integral terms describing the contributions of the 



1 In Ref. [5] we used the functions h? and Hj which differ from 
and H? by their normalization. The present form of the 
vertex functions is more convenient for the determination of 
x-dependent bare parameters. 



2 In Ref. p|] the quantity was denoted by J2. We changed 
the notation here in order to avoid its confusion with the 
notation of the electromagnetic current operator (see below). 
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three-body state to T 2 are 

ii(R±,x) = / R' ± dR' ± / dx' ^(-1/ 
Jo Jo . -,_ n 



t,i'=0 

(30a) 



Jj(12±,a;) = / flldfli / dx' ^ (-1/ 

Jo Jo { ,.,_ n 



t,i'=0 

c' ni h{(R' ± ,x') + C' ni Hf(R' ± ,x')_ , 

(30b) 



with n = 0, 1. The coefficients c, C, c ; , and C", which 
depend on R±, R' ± , x, x' ', j, and j 7 , are given in 
Ref. 0. 

We should now use the renormalization condi- 



tions rtl5b and (16) in order to determine the bare 



coupling constant g' 03 

terms of the functions h\ and H?, these conditions 
write 



(31a) 



h° (Rl,x) = g^l-g^\ 
HS(R* ± ,x) = gJl-g^\2-x), (31b) 



where R\ is defined by Eq. (Jl9|). S etting R± = R* ± 
and j = in Eqs. (|23a|) and fl23c[, we demand the 



relations (31) to be valid for arbitrary < x < 1. As 



explained in the previous section, this necessitates to 
substitute g' 03 — >■ ^03( x ) an d ^ — >• ^^(x). Using that 
at i?_L = i2* an d j = we have si = m 2 , and, hence, 
if? = z and K$ = 0, we get 



9o3K x ) = 



1 - 9 n\ 



2 r(2) 



-z°(i?l,x) 



(32a) 



9o3( x ) z L(x) = g 



/2 r(2-x)z[](i?l,x)-/ (i?l, a : 



1 -x 



(32b) 



Substituting these quantities back into the system of 
equations ( 23 ) , we obtain 



h J (R ± ,x) 


= V9 


+ 


,9' 3 


h{(R±,x) 






,9' 3 


H&(R±,x) 


= rjg(2-x) 


+ 


,9' 3 


H((R±,x) 


= V9 


+ 





Kihi(R ± ,x)+K J 2 h{(R ± ,x) 
-K 3 3 hi(R ± ,x) + K{h{{R^x 
K{Hi{R^x)+KiH{{R^x) 
-K J 3 H J (R ±1 x) + K{H{{R^x 



+ g^Ai^R^x), 
+ ^ /2 zj(i? ± ,x), 
+ g' 2 AI^(R±,x), 
+ ^ 2 Al{{R^x), 



(33) 



where 



and 



77 = 



l + <? 2 ^ 



(34) 



Ai 3 (R±,x) = z J (i? ± ,x)-z-g(i?l,x), (35a) 
A^'(i? ± ,x) = li(R ± ,x)-lS(R* ± ,x), (35b) 
Al((R ± ,x) = ll(R ± ,x)-i%(R* ± ,x). (35c) 

The functions /i] and Hf, being a solution of the in- 



homogeneous system of equations (33), are properly 
normalized. Once we have carried out the renormal- 



ization procedure [resulting in Eqs. (32)], these equa- 



tions contain only the physical coupling constant g, 
since g' is expressed via g by Eq. ( 25 ) . Let us remind 



that, in contrast to the results reported in Ref. [5], we 
do not fix here any particular value of x at which 
the renormalization conditions are considered. We 
keep the full ^-dependence of the bare coupling con- 
stant and the (^-dependent counterterm, according to 
Eqs. ([321). 



A Equations for the Fock components 



10 



Each index i and j can take two values, so that 
we have to deal with eight vertex functions. The 
convergence of the integrals over dR' ± in the integral 
terms (30) is ensured by the mutual cancellation of 



the physical and PV components at R' ± — > oo due to 
the following properties: 

h>i(R± ->• oo, a;) = hj(R± ->• oo,x), 



H^(R ± -> oo, x) = Hl(R± ->• oo,x), 



which automatically follow, when ^> fi\, from the 
structure of the equations ((331). This is a direct con- 



sequence of the PV regularization scheme. However, 
if one tries to go over to the limit ji\ — > oo inside the 
integrals (i. e. before the integration over dR' ± ), some 
of them become divergent, unless all the functions h\ 
and H 3 - vanish, when R' ± — > oo. But because of the 



nonzero free part on the right-hand side of Eqs. (33) 



there is no any reason to expect such a behavior of 
the solution. For instance, assuming that h\ and Hf 
at R' ± of order \i\ or greater have some finite val- 
ues, we would encounter divergences like log/ii. The 
renormalization leading to the subtractions (35) also 



does not fully prevent from divergences, because 
it "improves" the asymptotic (R' ± — > oo) behavior of 
the integrands only in the two terms, Ai^R^.x) and 
AIq(Rj_,x) (corresponding to the components with 
the physical external particles) among the eight ones. 
As a result, we can not take the limit /ii — >• oo directly 
in the system of equations p3| ), and the solution es- 
sentially depends on fii. 

From a practical point of view, we should solve the 
equations at large but finite values of \i\ , express phys- 
ical amplitudes (in which, by definition, all the exter- 
nal lines correspond to physical particles) through the 
vertex functions, and calculate observables. We then 
repeat these steps, gradually increasing If the cal- 
culated observables, as a function of /ii, tend to stable 
values within the required accuracy of the numerical 
calculations, we shall conclude that our renormaliza- 
tion scheme is successful. Let us emphasize that the 
stability is analyzed with respect to observable quan- 
tities only, while the vertex functions, as well as the 
Fock sector norms may strongly depend on fii , even if 
it reasonably exceeds the physical masses. In the fol- 
lowing, we will make this procedure numerically, for 
the calculation of the electromagnetic form factors. 



B. Calculation of the electromagnetic form 
factors 

The general decomposition of the spin- 1/2 electro- 
magnetic vertex (EMV) in CLFD is given by [3H51IT5] 



u(p')G p u(p) = eu(p' 

+ £; 



2m 

uj-p J uj-p 



' (uj-p) 2 _ 



u(p), 



(36) 



where P = p + p', q = p' -p, g? v = i(Yl v -l"l p )fa 
e is the physical charge, Fi and F 2 are the physical 
form factors, and £1,2,3 are nonphysical contributions. 
These latter originate from possible breaking of rota- 
tional symmetry, caused by the Fock space truncation. 
Under the condition u-q = 0, all Fi, 2 , £1,2,3 depend 
on Q 2 = — q 2 only. The physical form factors can be 
found according to [15] 



eF 1 
eF 2 



Tr + m)u p G p (tf + m)</i\ 



S(uj-p) 2 



(37a) 



m 



2{u-p)Q 2 



xTr 



(37b) 



Explicit analytical expressions for the two form fac- 
tors F\ and F 2 in the Yukawa model for the case of 
the / + fb + fbb Fock space truncation are given in 
Ref. [5]. For this reason, we will not dwell on techni- 
cal derivations, but focus on how the FSDR scheme 
works to renormalize the fermion EMV. 

From Eqs. (37) we can see that both physical 



form factors are determined by the contraction of the 
EMV with the four- vector u, sandwiched between the 
bispinors. It is thus convenient to define the operator 



eJ(Q) = 



u(p')u p GPu(p) 
2(u-p) 



(38) 



In standard LFD it is nothing else than the plus- 
component of the electromagnetic current. Within 
the FSDR scheme in truncated Fock space, the op- 



erator (38), similarly to the vertex functions, should 



be supplied with a superscript indicating the order of 
truncation N. J^ N \Q) can be represented as a super- 
position of contributions ji N \Q) from Fock sectors 
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with different numbers n of particles: 

N 



eo(N-n+i)jL N) {Q), (39) 

n=l 

where e (7v- n +i) is the electromagnetic bare coupling 
constant in the n-body sector [5 . By definition, 
eoi = e. The standard renormalization condition for 
the EMV 

G p \ p/=p = eY (40) 

written in terms of J(Q) has a very simple form 

J(0) = 1. (41) 

If the norms I n of all Fock sectors are finite, the fol- 
lowing relation must be valid: 



(42) 



for any N. The normalization condition ([8j just guar- 
antees the property d4lL In the Yukawa model how- 



ever, the norm of each Fock sector is infinite, and 
the validity of Eq. ( |42| ) depends on the regularization 
scheme. In Ref. [4J it was proved that Eq. (42) held 
true for any Fock sector containing one fermion and 
arbitrary number of bosons, provided the PV regular- 
ization is used, when mi — >• oo at finite fi\. Since this 
is just the case we consider in the present section, we 
can safely use the relation (42). A direct consequence 



of the latter is the following result [4]: 



eon = e, 



(43) 



for any n, i. e. the electromagnetic coupling constant 
is not renormalized at all. 



' e 03 



f\ 



FIG. 3: One-body (/) Fock sector contribution to the elec- 
tromagnetic vertex. 

In the / + fb + fbb approximation, the EMV is 
given by a sum of contributions shown graphically in 
Figs.[3}|5] The diagrams in Figs. [3] and [4] correspond, 
respectively, to the one- and two-body sector contri- 
butions to the EMV. The three-body sector contribu- 
tion in Fig. [5] is expressed, for convenience, through 
the two-body vertex and separated into two parts (a) 



' e 02 



FIG. 4: Two-body (fb) Fock sector contribution to the 
electromagnetic vertex. 



1 
1 

902 1 e i go2 




(a) 



(b) 



FIG. 5: Three-body (fbb) Fock sector contributions to the 
electromagnetic vertex, expressed through the two-body 
Fock component, with a nested fermion-boson loop (a) 
and with crossed boson loops (b). 



and (b). We will denote these particular contributions 
to the full operator J(Q) as J\(Q), J 2 (Q), J 3a (Q), and 
J?>b{Q)- To shorten notations, hereafter we omit the 
superscript TV for all quantities related to the three- 
body approximation, keeping it however, if a given 
quantity is calculated within a truncation of another 
order. From Eq. ( 39 ) we find 



eJ(Q) = e 3/i+e 2 J 2 (Q)+e [J 3a (Q) + Jsb{Q)] ■ (44) 



Due to Eqs. (41) and ( |43| ), the renormalized J(Q) be- 
comes 

J(Q) = 1 + [J 2 (Q)- J 2 (0)] 

+ [ha(Q) + J 3 b(Q) - J3a(0) - J 36 (0)] • 

(45) 

It is convenient to separate the common normaliza- 
tion factor \J 1 — g 2 I^ (I^ logarithmically diverges, 
when ji\ — y 00) from the vertex functions by defining 
a new set of functions 



1 



(2) 



(46) 



We will supply with a bar each contribution to J(Q), 
calculated through the functions and Hf in the 
vertices. In J 3a (Q) and J 3 b(Q) we suppose that the 
factor <5Tq 2 coming from the internal vertices in Fig. [5] 
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is also separated. Eq. (45) transforms thus to 



canceled in the difference 



J iff) = l + g 2 (l-g 2 li 2) )[J 2 (Q)-J2(0)] 



(47) 



x [J 3 „(Q) + J 3 t(Q) - J 3 a(0) - J 36 (0)] . 



Substituting here g 2 2 from Eq. ( 29a ) and rearranging 
the order of terms, we obtain 

Jiff) = i + g 2 [MQ)~ <Mo)] 



W{[j3a(Q)-I^J2(Q) 
J3a(0)-4 2) J 2 (0)]} 

W [JsbiQ) - J 36 (o)] • 



(48) 

Let us emphasize that although Eq. p8| ) looks like a 
perturbative expansion, it has no relation to perturba- 
tion theory, since the quantities J2, J^a, and J35 have 
rather complicated dependence on the coupling con- 
stant g, governed by the nonperturbative equations 
for the vertex functions. Nevertheless, it is rather 
instructive to study how fi\- divergences disappear in 
Eq. (48), when the vertex functions are given by their 



perturbative values: 

h{ = H{ = 1, h{ =0, H 3 =2- x. 

In this case, the fi\ -dependence of J2, J^a, and J36 
comes only from the integrals for the EMV. It is easy 
to show that both J 2 and J-^b diverge, when \i\ — >> 00, 
as log //1 with a coefficient independent of Q. Sub- 
tracting from each of them its value at Q = 0, we just 
cancel divergent logarithmic terms, and the result is 
finite. Concerning the contribution J 3a [see Fig.jHJa)], 
it has a nested fermion-boson loop which diverges, af- 
ter the integration over the corresponding kinematical 
variables, as log The integration over the variables 
of the external loop gives one more log/xi, so that 
J 3a diverges like log 2 pL\ . Its renormalization occurs 
in two steps, as is seen from the expression in braces 
in Eq. ( [48| . In a first step, we form the difference 
— —(2) ~ 1 

Jd,a(Q) — I2 JziQ) which kills log/ii coming from 

the nested loop. Indeed, by definition, I 2 coincides, 
up to a factor of g 2 , with the two-body normalization 
integral which, according to Eq. ( [42] ) , is related to the 
two-body contribution to the EMV by 



(o) = 4 2) . 



(49) 



Due to the fact that the divergent part of the nested 

—(2) 

loop coincides with that in J 2 (0), it is completely 



J3a(Q)-I^HQ) 



In a 



second step, subtracting from this difference its value 
at Q = 0, we remove divergences coming from the in- 
tegration over the external loop variables. In standard 
perturbation theory, the renormalization is made by 
the same scenario. But we emphasize that Eq. (48) 
naturally appears within the FSDR scheme which is 
fully nonperturbative. 

When the vertex functions are solutions of the sys- 
tem of equations ( 33 ) and depend non-trivially on par- 



ticle momenta and on the coupling constant, it is much 
more difficult to trace analytically the cancellation of 
Hi- divergences in calculated observables. For this rea- 
son, we will make the corresponding analysis numeri- 
cally. 



C. Numerical results 



We solve the system of integral equations (33), 



in order to calculate the two-body Fock component. 
Knowing the latter, we are able to find the three-body 
component as well (see Ref. [5] for details). After 
that, we calculate the electromagnetic form factors 
by means of Eqs. ([37]), with the EMV G p given by a 
sum of contributions shown in Figs. |3H5l The EMV 



is expressed through the quantity J(Q) for which we 
take its renormalized value ( |48) . In all computa- 
tions, we use the physical particle masses m = 0.938 
and fi = 0.138 reflecting the characteristic nuclear 
physics mass scales. Each physical quantity is cal- 
culated for three values of the physical coupling con- 
stant, a = g 2 /An = 0.5, 0.8, and 1.0. 

We first compute the fermion AMM which is the 
value F2(Q 2 = 0). It is shown in Fig. [6] as a function 
of the PV boson mass fi\. Each of the two- and three- 
body Fock sector contributions to the AMM essen- 
tially depends on while their sum is stable, as \i\ 
becomes large enough. Note that using x-dependent 
bare parameters removes /ii-dependence of the AMM, 
observed in Ref. [5] already for a ~ 0.5, even for larger 
coupling constants. 

In Figs. [7] and [8[ we plot the electromagnetic form 
factors F\ and F 2 , respectively, as a function of Q 2 for 
Hi = 100, for the three values of the coupling constant, 
considered here. For this value of /ii, the form factors 
already reach the zone of stability: further increase 
of fii does not lead to changes distinguishable in the 
scale of the figures. 

We finally show in Fig. [9] the contributions of the 
one-, two-, and three-body sectors to the norm of the 



C Numerical results 



13 



0,125 
0,100 



^ 0,075 
O 



0,050 
0,025 - 
0,000 



a = 0.5 



5 6 7 

log 10 l>i 2 /n 2 ] 



0,20 



0,15 



r> °' 10 

o 



0,05 
0,00 

0,25 

0,20 

^ 0,15 
O 

(N 

U. 0,10 
0,05 
0,00 



a=0.8 



5 6 7 

log 10 l>i 2 /n 2 ] 



a=1.0 



3 4 5 6 7 

logio [niW] 

FIG. 6: The anomalous magnetic moment in the Yukawa 
model as a function of the PV mass fii, for three differ- 
ent values of the coupling constant, a — 0.5 (upper plot), 
0.8 (middle plot) and a — 1.0 (lower plot). The dashed 
and long-dashed lines are, respectively, the two- and three- 
body contributions, while the solid line is the total result. 
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FIG. 7: Electromagnetic form factor F\ (Q 2 ) in the Yukawa 
model, at fii — 100, for a — 0.5 (upper plot), 0.8 (mid- 
dle plot) and a — 1.0 (lower plot). The dotted, dashed, 
and long-dashed lines are, respectively, the one-, two-, and 
three-body contributions, while the solid line is the total 
result. 



state vector. It is interesting to note that for large 
enough values of the PV boson mass fi\, the norm 
of the one-body Fock sector gets negative, while the 
norm of the two-body Fock sector becomes larger than 
1, the sum of all contributions being fixed to 1 by 
construction. However, physical observables are well 
denned and do not show any discontinuity, as a func- 
tion of Hi, when Ii changes the sign or I 2 exceeds 
unity. This is an illustration of the fact that the 
norms of the Fock components, in the presence of the 
PV sectors having negative norms, are not physical 
observables and, hence, they are expected to be reg- 



ularization scale dependent. As we can see in these 
figures, they do depend on fii, unlike the electromag- 
netic form factors. The stability of the form factors to 
variations of the PV masses, when the latter ones are 
large compared to the physical masses, opens encour- 
aging perspectives for our method of nonperturbative 
renormalization. 

For obvious reasons related to the finite accuracy 
of numerical calculations, this domain of stability 
can not be checked up to infinite values of the 
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FIG. 8: The same as in Fig. [7J but for the ratio 
F 2 (Q 2 )/F 2 (0). 
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FIG. 9: Individual contributions of the one- (dotted line), 
two- (dashed line), and three-body (long-dashed line) Fock 
sectors to the norm (solid line) of the state vector, as a 
function of the PV boson mass /xi, for a — 0.5 (upper 
plot), 0.8 (middle plot), and a = 1.0 (lower plot). 



PV mass fi\. However, we have found that for 
enough large fi\ the stability disappears. Namely, 
when fjb\ approaches some critical mass \i\ Cl all the 
functions h\ and Hj defining the two-body Fock 
component become very large, when R± — >• oo or 
x — >• 0. At fi\ = (i\ c these functions are unbounded, 
because their asymptotic values turn into infinity. 
Since hj, Hj start growing from the characteristic 
values R± ~ fi\ and x ~ (m//ii) 2 , the calculated 
observables are sensitive to the variations of /ii, if 
it is not far enough from \i\ c . The critical mass 
/iic which can be calculated analytically has very 
sharp dependence on the coupling constant a. For 
a — >■ it goes to infinity. As a increases, fii c falls 



down very rapidly. Thus, for the physical masses 
we consider (m = 0.938, fi = 0.138) the values 
of fiic, corresponding to a = 0.5, 0.8, and 1, are, 
respectively, 5.17-10 8 , 4.17-10 4 , and 1.8-10 3 . For little 
bit larger en's, fii c becomes comparable with the 
physical masses. In this domain, a reliable numerical 
calculation of physical observables is impossible. 
So, in the Yukawa model, within the three-body 
truncation and using the PV regularization scheme, 
the typical dependence of calculated observables on 
fii looks as follows. At fii comparable with the 
physical masses, observables are sensitive to the 
value of fii. Then, in the region m <^C fii <^C fiic 
this dependence has a plateau corresponding to the 
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zone of stability. This is just our "working" region, 
inside which we can trust the physical meaning of 
our numerical results. At (i\ close to fi\ c we can 
not expect to perform reliable numerical calculation 
of physical observables. The determination of the 
critical mass /ii c , as well as studying the properties of 
the system of equations ( p3| , when \i\ is close to \i\ Cl 
go beyond the framework of the present article and 
will be explained in detail in a future publication. 



from Eqs. (|33). 



D. 



Additional test of the renormalization 
scheme 



The system of equations (33) was obtained in the 
limit mi — > oo, while \i\ being fixed. We then in- 
crease fii till the stability of calculated observables is 
reached. Physical observables however must be inde- 
pendent of the order in which the infinite PV mass 
limit is taken. For this reason, we shall make a test of 
the self-consistency of our renormalization procedure. 
Namely, we perform an additional calculation of the 
AMM, taking first the limit fii — > oo (analytically) 
and then the limit mi — > oo (numerically). 

The corresponding system of equations for the ver- 
tex functions in the limit fJL\ — > oo can be obtained, 
analogously to Eqs. p3| ), from the general equation 
shown graphically in Fig. [2] Omitting here all tech- 
nical details, we indicate only the main differences of 
the system of equations obtained in the limit \i\ — >• oo 



When fii — > oo, the four vertex functions with j = 
form a closed sub-system of linear integral equations, 
while the other four equations involving the functions 
with j = 1 can be omitted, since the latter ones do 
not contribute to the AMM. 

As far as the calculation of the electromagnetic form 
factors is concerned, the change of the order of PV 
mass limits also brings some new features to the pro- 
cedure. Eq. (42), and following from it Eq. (p3|, are 
not anymore valid. The same relates to Eq. (49|) which 
is a particular case of Eq. p2| >. As a result, we can 



not rely on the renormalized expression (48). Indeed, 



(2) 

the norm of the two-body sector m the two-body 
approximation was calculated in Ref. [4] [see Eq. (49) 
there]. Amputating from it the factor g 2 , we find 



f(2) _ 
±o — 



1 

8^ 



/ dR±R± dxx V(-l)* +J 
Jo Jo ij=0 

R\ + [(1 -x)m + mi] 2 
[R 2 ± + (1 - x)\i 2 + xm 2 - x(l - x)m 2 } 2 ' 

(50) 



—(2) 

Whereas, for the electromagnetic current J2 (0) = 
J2 2 \o)/g 2 , where is the two-body component of 
defined by Eq. H with the EMV G p found in the 



two-body approximation, the calculation gives (the 
details can be found in Ref. [5]): 



Jf(0) = 



1 

8^2 



/ dR± R± dxx Y] (-1) 
./o Jo 



i+i'+j 



,j=0 

R\ + (1 — x) 2 m 2 + (1 — x)m(rrii + m^) + 



[R 2 ± + (1 - x)ii 2 + xmf - x{l - x)m 2 ][R\ + (1 - x)/2 2 + xm 2 , - x{l - x)m 2 ] ' 



(51) 



Both integrals (50) (the two-body norm) and (51) limit the difference j!f\o) —ffl is not zero. For large 



(the two-body current at zero momentum transfer) 
converge in the limits mi — > 00 for finite fii and 
lii — >• 00 for finite mi. As mentioned above, in the 
limit mi — > 00 for finite /ii, they coincide with each 
other and satisfy the relation ( |49| . Taking the op- 
posite limit yL*i — y 00 for finite mi means retaining 
only the terms with j = (the physical boson index) 
in the sums over j. One can easily check that in this 



values of mi, it tends to a finite, mass independent, 
value: 



(2) (0)-7f = 



167T 



+ 



m n mi 
— log — 
mi m 



(52) 



while both J2 (0) and 7 2 taken separately diverge 
as log(mi/m), when mi — > 00. So, under the PV reg- 
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ularization scheme, the contribution of a particular 
Fock sector to the charge form factor at zero momen- 
tum transfer may not coincide with the norm of this 
sector. 

To get a renormalized expression for J(Q), we re- 



peat the same steps as in Sec. |IVB] but without using 
Eqs. (43) and (49). From the renormalization condi- 
tion 



( |4ip considered, by turn, in the two- and three- 
body truncated Fock spaces, we determine the bare 
electromagnetic coupling constants eo2 and eo3- Now 
both of them differ from the physical charge e. Sub- 
stituting them into the general formula ( |44| ), we arrive 
at the result for the renormalized J(Q), which differs 



from Eq. (|48f by the substitution J^ 2) J^ Zj (0). 



f( 2 )/ 



We compare in Table |T| the numerical results for 
the AMM, obtained in this way, with those found in 
Sec. [IVC] in the limit mi oo. The AMM is consid- 
ered as a function of the PV mass which is kept finite 
(mi, if the limit fi\ — > oo has been taken, and vice 
versa). For convenience of the comparison, we took 
the same sets of finite PV mass values. 



V. ANTIPARTICLE DEGREES OF 
FREEDOM 

A. Contribution to the two-body vertex 
function 

We can extend the Fock decomposition of the 
fermion state vector by introducing the antifermion 
d.o.f. In the lowest (also three-body) approximation 
this corresponds to adding the /// Fock sector to 
those previously introduced (/, fb, and fbb). We have 
already considered the role of the three-body Fock sec- 
tor with an antiparticle within the pure scalar model 
(a heavy scalar boson interacting with light scalar 
bosons) in Ref. [6]. We perform here a similar study 
in the Yukawa model. 

The antifermion d.o.f. contributions to the two- 
body vertex are of the following two types: 

(i) The first one corresponds to standard fermion- 
antifermion polarization corrections to a boson line, 



as shown in Fig. 10 



TABLE I: The anomalous magnetic moment calculated 
for a — 0.8 in the two different limits of the PV masses. 



PV mass kept 


AMM when 


AMM when 


finite (fii or mi) 


mi — > oo 


fll — > oo 


5 


0.1549 


0.1454 


10 


0.1641 


0.1630 


25 


0.1690 


0.1704 


50 


0.1702 


0.1715 


100 


0.1706 


0.1716 


250 


0.1708 


0.1714 


500 


0.1709 


0.1713 



If each of the finite PV masses is much larger than 
all physical masses, the values of the AMM, obtained 
in both limits, coincide within the computational ac- 
curacy (about 0.2%), as it should be if the renormal- 
ization procedure works properly. We can thus choose 
any convenient order of the infinite PV mass limits. 
Since the equations for the Fock components are tech- 
nically simpler in the limit mi — >• oo, we continue 
working with the vertex functions and the EMV taken 
in this limit. The independence of physical results on 
the order in which the infinite PV mass limit is taken 
and, hence, on the way we use to get rid of the bare 
parameters, is a strong evidence of the self-consistency 
of our renormalization scheme. 




FIG. 10: Contribution of the /// Fock sector to the two- 
body vertex: polarization correction to the boson line. 
The antifermion is shown by the thick line. 

(ii) The second one corresponds to transition am- 
plitudes fb — > fb associated to the excitation of an- 
tifermion d.o.f. from a fermion line, as shown in 
Fig. Ha). 





(b) 



FIG. 11: Contributions of the /// (a) and fbb (b) Fock 
sectors to the two-body vertex due to the transition am- 
plitude fb — > fb. 

These two contributions (i) and (ii) have different 
nature. The first one is the // loop. The second con- 
tribution is the fermion-boson loop. In addition, the 
diagram in Fig. [rT|a), by changing the order of ver- 
tices with respect to the light-front time, evolves to the 
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fbb Fock sector contribution, as shown in Fig. [TTj^b). 
We shall consider, in this paper, the contributions 
of the second kind only, which corresponds to the 
quenched approximation in LFD. 

The system of equations for the Fock components 
can be obtained by direct generalization of the pro- 
cedure exposed in Ref. [6j for the scalar case. We in- 
troduce one more three-body (///) Fock component, 
in addition to the /, fb, and fbb ones considered in 
Sec. 



IV 



In the three-body approximation, this new 
Fock component is easily expressed through the two- 
body component, as well as the fbb one. As a result, 
we obtain a closed (matrix) equation for the two-body 
vertex function, as given, in the quenched approxima- 
tion, by Fig. [12] It differs from the equation in the 
/ ' + fb + fbb approximation, shown in Fig. [2] by an ad- 
ditional term on the right-hand side (the last diagram 



in Fig. 12). 



We should pay attention to the fact of using the 
same constant, #02, in the elementary vertices / ^ 
/ + 6 and b <r+ f + f. Strictly speaking, this is not 
mandatory from the point of view of the general FSDR 
rules, because contributions from fbb and /// states 
represent different Fock sectors. Nevertheless, we as- 
sign to these vertices the same factor g 02 , which seems 
quite natural in this first study of the influence of an- 
tifermion d.o.f. 



since we do not consider in this study polarization 
corrections to boson lines, the renormalization condi- 
tion Jl6b is untouched. The second condition Jl5j also 



does not change, because it is universal. We can thus 



proceed further in the same way as in Sec. IV A The 
form of the system of equations for the functions h 3 
and H 3 remains the same as in Eqs. (33), but the in- 



tegral terms i 3 and I 3 Q [see Eqs. (30)] obtain additive 
j-independent contributions according to 



i 3 ^i 3 + i , I J ^I 3 +I , 



(53) 



where 
io(R±,x) 



/ R' ± dR' ± / dx' V(-l/ 

JO Jl-x j/=0 

x 'c Q h{{R'^x') + C Q H 3 Q \R'^x')] , 



(54a) 



I (R±,x) = 



/ R' ± dR' ± / dx' 

JO Jl-x j/=0 

x ~c' h((R' ± ,x') + C' Hi'{R' ± ,x') . 



7(3) r 



(54b) 

The integral terms i{ and I[ do not change. The co- 
efficients Co, Co, c' , and C' Q determining antifermion 
contributions are given in Appendix [Aj Note that the 
limits of the integration over dx' in Eqs. (54) differ 
from those in Eqs. (|30). 






FIG. 12: Graphical representation of the equation for the 

(3) 

two-body vertex function including the contribution 
of antifermion d.o.f. in the quenched approximation. 



Since the general structure (12) of the two-body 



vertex function in the Yukawa model is universal, and 



B. Numerical results 

We have solved numerically the system of equa- 
tions (33) with the integral terms modified accord- 



ing to Eq. (53), for the same set of parameters as in 
Sec. |TVC] m = 0.938, \i = 0.138, and a = 0.5, 0.8, 
and 1.0. Along with the functions h 3 and Hf, we cal- 
culate al so th e bar e par ameters g' 03 and Z' u defined 
by Eqs. (26a) and (26b), respectively. In truncated 



Fock space, both of them are functions of x, accord- 



ing to Eqs. (32), where the integral terms include now 



antifermion contributions. Besides that, they depend 
also on the PV mass 



We plot in Figs. [13] and [14] these bare parameters 
as a function of x, each for a = 0.5, 0.8, and 1.0, at 



a typical value fi\ = 100. In Fig. [13] the relative value 
of #03 with respect to its mean value g' ' 03 over the 
interval < x < 1 is shown, i.e. we plot the quantity 



Wo3( x ) = Wosix) ~ £03] /to 
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where g' Q3 = J* g' Q3 (x)dx. For comparison, we show 
also on these plots the same functions calculated with- 
out antifermion contributions. The most interesting 
fact is that the function g' 03 (x), which exhibits strong 
^-dependence in the / + fb + fbb approximation, be- 
comes almost a constant, if the /// Fock sector is 
included. Concerning the function Z' u (x), it shows 
a similar tendency as well, with a bit stronger x- 
dependency than g' 03 (x). In addition, the magnitude 
of Z' u (x) is reasonably smaller than that calculated in 
the / + fb + fbb truncated Fock space. 



0,4 



0,2 - 



0,0 - 



-0,2 



a=0.5 



0,0 0,2 0,4 0,6 0,8 1,0 



0,4 



0,2 



0,0 



-0,2 



a=0.8 



0,0 0,2 0,4 0,6 



0,8 



1,0 



0,4 



0,2 



0,0 



-0,2 



a=1.0 



0,0 0,2 0,4 0,6 0,8 1,0 



FIG. 13: x-dependence of the bare coupling constant g' 03 , 
calculated relatively to its mean value over the interval 
x £ [0, 1], for a — 0.5 (upper plot), a — 0.8 (middle plot) 
and a — 1.0 (lower plot), calculated for fn — 100. The 
solid (dashed) lines correspond to the results obtained with 
(without) the /// Fock sector contribution. 



a=0.5 



0,0 0,2 0,4 0,6 0,8 1,0 



6 - 
4 - 



-f 



a=0.8 



0,0 0,2 0,4 0,6 0,8 1,0 




0,0 0,2 0,4 0,6 0,8 1,0 



FIG. 14: x-dependence of the counterterm Z'^ for a — 0.5 
(upper plot), a — 0.8 (middle plot), and a — 1.0 (lower 
plot), calculated for \i\ = 100. The solid (dashed) lines 
correspond to the results obtained with (without) the /// 
Fock sector contribution. 



stants is not a specific property of the Yukawa model, 
since we already encountered similar features in our 
studies of the scalar model [6^ . In contrast to the lat- 
ter, where we incorporated particle-antiparticle loop 
contributions as well, we reproduced the property 
#03 ( x ) ~ const in the Yukawa model within the 
quenched approximation. From here it follows that 



namely the contribution shown in Fig. 11 a) is respon- 



sible for this property. This was the main reason of us- 
ing here the quenched approximation which allows us 
to keep such an important property of the bare param- 
eters, on the one hand, and to retain the renormaliza- 



The fact that gb 3 (x) and Z' u (x) are close to con- tion condition ( 16 ) and avoid complications connected 
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with the fermion-antifermion loop renormalization, on 
the other hand. 

For a deeper understanding of what happens with 
the bare parameters, when the antifermion d.o.f. are 
involved, let us consider the following simple approx- 
imation. It is instructive to solve the equation in 



Fig. [12] by iterations, taking, as the zero order ap- 
proximation, the value of the two-body vertex func- 
tion obtained for the f + fb Fock space truncation, 
i. e. with T 2 = g. We then substitute this value on the 
right-hand side of this equation and calculate the first 
iteration, T 2 (R±,x), in the limit mi — > oo. We can 
then calculate T 2 (Rj_(x),x) [see Eq. (fig])], with both 



constituent particle legs corresponding to the physical 
particles (i = j = 0). This just determines the func- 
tions bi j2 (R]_(x) 1 x) which enter into the renormaliza- 
tion conditions ( [15] ) and (16). Since we are interested 
in the x-dependence of these functions, it is enough 
to calculate the contributions to the latter ones from 
the diagrams shown in Fig. [TT] By direct calculation, 
we find that the sum of these two contributions to the 
function b\ (denoted by bf and b\) on the energy shell 
does not depend on x: 

b%(R* ± (x),x) + b b 1 (R* ± (x),x) = const, (55) 

while the total value of b 2 is zero: 

bS(I? ± (x),x) + b*(R* ± (x),x) = 0. (56) 

So, on the level of the first iteration, we can meet 



the renormalization condition (16) exactly, with g' 03 
independent of x. The other condition (PL5|) is satis- 



fied automatically, without the need of an additional 
counterterm (20), i.e. Z u (x) = 0. Note that the coun- 



terterm defined by Eq. (26b) does not turn into 
zero, because it includes, beside Z u , the contribution 
from the one-body state. We can assert only that in 
the approximation discussed above it becomes a con- 
stant. 

This result has a simple explanation. The diagrams 
shown in Fig. [TT] with a constant internal two-body 
vertex T 2 = g coincide with the light-front perturba- 
tive ones taken in the order g 3 . Since there are no 
other g 3 -order perturbative contributions to the two- 
body vertex function, the sum of their amplitudes, 
on the energy shell, is identical to the corresponding 
on-mass-shell Feynman amplitude which is a constant 
and does not depend on u. Eqs. (55) and (56) are 
direct consequences of this fact. 



We may continue iterating the equation in Fig. [12 
and represent each of the functions b\^ 2 as a series 



in powers of the coupling constant. But these ex- 
pansions, starting with the order g 5 , differ from the 
perturbative ones. Indeed, perturbative contributions 
involve Fock sectors with arbitrary number of parti- 
cles (of course, only those which are compatible with 
the order of perturbation), while the nonperturba- 
tive equation for T 2 contains no contributions from 
higher than three-body Fock sectors. As a result, 
^i,2(^l(^), x) are no more constants, i. e. they de- 
pend on x. In order to enforce the fulfilment of the 
renormalization conditions, we have to introduce x- 
dependent bare parameters g' 03 (x) and Z'^{x). The 
same happens in our nonperturbative calculations, 
where the two-body vertex function T 2 is far from 
being a constant. Note that if we assigned different 
vertex factors to the elementary vertices / — »• / + b 
and b — >• / + /, we would not get the properties (55) 
and J56). 



In the f+fb+fbb approximation, the ^-dependence 



of b\ 



_{x),x) 



is governed by the diagram in 
Fig. [lllb) only. The lowest order iterative contri- 
bution of the latter diagram to the function b 2l at 
/ii ^> {m, /i}, has the form 

bURl (x),x) = -£ 2 - log ^ + f 2 (x) , (57) 



where f 2 (x) is a function of x and it does not depend 
on In order to save the renormalization condi- 



tion (15), we have to add the structure (20) to the 



interaction Hamiltonian, with the ^-dependent coun- 
terterm Z UJ (x) = —b 2 (R* ± (x) , x) which has a con- 
stant /independent part (divergent when [i\ — > 00) 
and a finite ^-dependent part. For higher order it- 
erations, terms divergent when fi\ — > 00 appear in 
the x-dependent part of the counterterm as well. 
Hence, in our nonperturbative calculations within the 
/ + fb + fbb truncated space, we do expect strong 
x-dependence of the bare parameters, which is con- 



firmed by the results represented in Figs. 13 and 14 



We can now proceed to the calculation of the elec- 
tromagnetic form factors. The general method is the 
same as in Sec. IIVB1 With the antifermion d.o.f. in- 
cluded in the quenched approximation, one should 
take into account the additional contributions to the 
EMV, shown in Fig. [15] The amount of computations 
can be reduced by exploiting some symmetry prop- 
erties of the diagrams, since the two contributions in 
Fig. [l5|a) are exactly the same. So, we can calculate 
only one of them and multiply the result by a factor 
of 2. 
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(«) 




FIG. 15: Antifermion contributions to the electromagnetic 
vertex in the quenched approximation. 



For numerical calculations, we take the same values 
of the particle masses and the coupling constant as in 
Sec. lIV Cl Our results for the fermion AMM are shown 



in Fig. 16 Although the antifermion d.o.f. have dras- 
tic influence on the ^-dependence of the bare param- 
eters, they decrease the AMM value, as compared to 
that in the / + fb + fbb approximation (see Fig. [6]) , 
by only ~ 2% for a = 0.5, ~ 4% for a = 0.8, and 
~ 8% for a = 1.0. So small influence of the /// Fock 
sector on the AMM is caused by the smallness of the 
ratio fi/m = 0.147. At larger fi the effect of the an- 
tifermion d.o.f. reveals itself stronger. The AMM is 
stable, as a function of /xi, if the latter is much larger 
than the physical masses, but not too close to the crit- 
ical mass iii c which is t he sam e as in the / + fb + fbb 
approximatoin (see Sec. IV C). If the latter condition 



is violated, we recover the AMM instability similar to 
the one obtained in the / + fb + fbb truncation, with 
however a larger amplitude. A visible deviation of the 
AMM from the constant value is seen on the lower 
plot in Fig. [THJ when fii > 200, while \i\ c for a = 1 
is about 1800 [log 10 (/if c //i 2 ) « 8.2]. For a = 0.5 and 
a = 0.8 [log 10 (fil c / fi 2 ) w 19.1 and « 11.0, respec- 
tively] the AMM deviations from the constant in the 
interval 100 < fi\ < 1000 do not exceed the computa- 
tional precision level (~ 0.5%). This means that with 
the same accuracy we have no any sign of uncanceled 
divergences. 
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FIG. 16: The same as in Fig. [6| but including the contri- 
butions from antifermion d.o.f. (dash-dotted line), shown 
in Fig. fT5l 



VI. CONCLUSION 

The results reported in this study are a first exam- 
ple of a full, nonperturbative, calculation of the prop- 
erties of relativistic compound systems in the FSDR 
framework. The general approach is based on an ex- 
pansion of the state vector of the system considered in 
Fock components, within LFD. The use of CLFD, to- 
gether with an appropriate renormalization scheme in 
truncated Fock space, gives a very promising oppor- 
tunity to calculate properties of compound systems in 
a regularization scale invariant way. 

The full implementation of the renormalization con- 
ditions, which relate the on-energy- shell two-body ver- 
tex function to physical observables, is the last build- 
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ing block in our strategy to get reliable predictions for 
physical observables in the nonperturbative domain. 
It leads a priori, and unambiguously, to the depen- 
dence of the bare parameters on one of the kinemati- 
cal variables of the two-body vertex function, like for 
instance the longitudinal momentum fraction x. This 
dependence must disappear in an exact calculation, 
i.e. when the Fock space is not truncated. 

We applied our strategy to the calculation of the 
fermion electromagnetic form factors in the Yukawa 
model, in the three-body Fock space truncation. 
Within our numerical precision the form factors, as 
a function of the momentum transfer, are indepen- 
dent of the regularization scale (the P V boson mass) , 
as soon as the latter is large enough compared to the 
typical intrinsic energy /momentum scales of the sys- 
tem, but smaller than some critical value (i\ c . 

We finally investigated the role of antifermion d.o.f. 
We showed that, in the leading order of perturba- 
tion theory, the contributions to the state vector of 
the Fock sector with an antifermion are precisely the 
ones which make the renormalization conditions fully 
consistent. In that case, no extra cu-dependent coun- 
terterms are needed to restore the rotational invari- 
ance of the two-body vertex function on the energy 
shell, while the x-dependence of the bare parameters 
is canceled exactly. In our nonperturbative calcula- 
tions, these contributions considerably improve the 
self-consistency of the renormalization conditions. 



and 



r 2n 



Jo = 



Ji = 



d<f>' 



sign(A) 



/o 2tiD(A + Bcos4)') DVA 2 - B 2 ' 
f 27r cosW 1 A - 



/o 2TrD(A + Bcosd>') DB B 
Now the coefficients in Eqs. (|54|) take the form 



c = 



Co = 



Cn 



./• - Dm. 



[R±R' ± Jo 
2x + 2x' - 

R ± (3x' - 2) J - R' ± (3x - 2) Ji] , 



R±m 

+m 2 (xx' + 2x + 2x' - 4) Ji] , 
(x - l)m 2 

R±m 



(x - 1)R' 



m 

(x-1) 



-\-R±R'j_Jlj . 



x [R' ± (3x - 2) J - R±(3x' - 2) Ji] , 



m 2 (xx / + 2x + 2x' -4) J 



Note that these coefficients do not depend on the in- 
dex j. 



Appendix A: The coefficients Co, Co, c'o, and C'o 

It is convenient to introduce the following notations: 

771 = R' 2 ± + x' 2 m 2 + {l-x')$,, 

A = (l-x)xR'\ + (l-x')x'R\ 

+xx'(2 — x — x')m 2 , 
B = 2(1 - x)(l - x')R' ± R±, 
D = 8tt 2 , 
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